Observation of Anomalous Phonon Softening in Bilayer Graphene 
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The interaction of electron-hole pairs with lattice vibrations exhibits a wealth of intriguing physical 
phenomena. The Kohn anomaly is a renowned example where electron-phonon coupling leads to 
non-analytic phonon dispersion at specific momentum nesting the Fermi surface 11] . Here we report 
evidence of another type of phonon anomaly discovered by low temperature Raman spectroscopy in 
bilayer graphene where the charge density is modulated by the electric field effect. This anomaly, 
arising from charge-tunable modulations of particle-hole pairs that are resonantly coupled to lattice 
vibrations, is predicted to exhibit a logarithmic divergence in the long-wavelength optical-phonon 
energy. In a non-uniform bilayer of graphene, the logarithmic divergence is abated by charge density 
inhomogeneity leaving as a vestige an anomalous phonon softening. The observed softening marks 
the first confirmation of the phonon anomaly as a key signature of the resonant deformation-potential 
electron-phonon coupling. The high sensitivity of the phonon softening to charge density non- 
uniformity creates significant venues to explore the interplay between fundamental interactions and 
disorder in the atomic layers. 
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Recently, the abihty to use the electric field effect 
(EFE) to continuously dope large densities of electrons 
or holes into graphene has led to general interest in 
the deformation potential coupling of optical phonons 
with charge carriers in the extreme two-dimensional(2D) 
atomic limit [1, 0, S H, 01 1 and to experimental obser- 
vations of charge-tunable phonon energy and lifetime in 
single layer graphene [1, [^, [lO| . In prior EFE-Raman ex- 
periments, however, a logarithmically-divergent phonon- 
energy anomaly Q predicted by the dynamical perturba- 
tion theory was not observed. The puzzling absence of 
the phonon anomaly was attributed to the presence of 
large non-uniformity of the charge density in monolayer 
graphene samples [8| . This interpretation is subsequently 
supported by the observation of asymmetric phonon line- 
shape in the EFE-Raman spectra at low doping levels fll'|. 

In this paper we demonstrate that the phonon 
anomaly, while proposed initially in single layer graphene 
0, is a generic feature of tunable resonant electron- 
phonon coupling when the particle-hole pair energy is 
close to the phonon energy (Fig.l). We report experi- 
mental observation of the anomalous phonon softening, 
a clear signature of such phonon anomaly, in bilayer 
graphene. The observation of the resonant phonon soft- 
ening is facilitated in bilayer graphene. The reason is that 



the coupling between the two graphene layers results in a 
nearly parabolic dispersion and relatively large density of 
states near the vanishing band gap, making the phonon 
anomaly robust even in the presence of relatively large 
charge density non-uniformity. 

Rather than being a special property unique to the 
monolayer graphene, the anomaly is anticipated to oc- 
cur in any multi-band electron system where the zero 
momentum optical phonon can create resonant electron- 
hole pairs across the conduction and valence bands. Fig- 
ure 1 illustrates the concept of the phonon anomaly in 
a 2D gapless semiconductor with particle-hole symmetry 
as a model system. A system with similar low energy 
band structure is in fact bilayer graphene, in which the 
two valleys of conduction and valence bands are centered 
at K and K' corner points of the Brillouin zone [l^ . 
The deformation-potential interaction couples vertical 
electron-hole interband excitations with long-wavelength 
optical phonons in graphene 1^ 14 1 . This coupling con- 
tributes to the renormalization of the phonon energy [l5| . 
Considering the electron-phonon interaction within 2nd 
order time-dependent perturbation theory [lB| , combined 
with the Pauli exclusion principle, we obtain the change 
of phonon energy Wp/jwith tuning of the Fermi energy 
Ep: 
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where Ee-h^s the energy of an electron-hole pair, A 
is the electron-phonon coupling parameter with dimen- 
sion of energy. Because of the resonant denominator 



ru,j h'^-E — ^ ™ equation (1), the perturbative contribu- 
tions to cjp/jhave a sign change at i?e_^=^a;p;jand are 
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FIG. 1: Illustration of the phonon anomaly concept, a-c, Vertical interband electron-hole pair transitions in a gapless 2D 
semiconductor with 3 different Fermi levels. Regions with green shading are filled with electrons. The transition indicated by 
the blue arrow is the resonance with the long-wavelength optical phonon. d, Predicted change of phonon energy as a function 
of the Fermi energy. The two phonon anomalies show up at Ef = :ttujjph/'2- 



markedly enhanced near the resonance (the interband 
transition indicated with a blue arrow in Fig. la). When 
\Ep\ is increased, the positive and negative enhanced 
perturbations to ujph are switched off due to restrictions 
placed by the Pauli Principle (illustrations for the case of 
electron doping are shown in Fig. lb,c) . For this reason, 
the modulation of carrier density in the system results 
in marked changes in the optical phonon energy around 
\Ef\ = ?ia;p/i/2, as shown in Fig. Id. 

Note that in the derivation of equation (1), no assump- 
tions of system dimensionality or band dispersion curva- 
ture were made. We expect the phonon anomaly in prin- 
ciple to show up in ID system like metallic or near zero 
gap semiconducting carbon nanotubes 13, graphene 
nanoribbons [3], 2D system like graphene thin films, 3D 
system like silicon and germanium [19, 2Q]- In fact, like 
in single layer graphene, similar log-divergence was pre- 
dicted, albeit not observed, in the studies of chemically 
doped silicon [19]. 

In the present work, we use EFE to tune Ep to access 
the resonant electron-phonon interaction condition of the 
phonon anomaly in bilayer graphene. Low temperature 
Raman spectroscopy was employed to detect the evolu- 
tion of the zero momentum optical phonon, known as 
G band in graphite and graphene thin films. Its energy 
hoja^ 1580 cm^^ w 196 meV. Bilayer graphene is an 
intriguing electronic system in which transport measure- 
ments have revealed unconventional quantum Hall effect 
and Berry's phase of 27r 2]J. It has 4 carbon atoms in 
its unit cell. The locations of them are At Bt in the top 
layer, Ai, B\j in the bottom layer (Fig. 2a). The pz orbital 
of these atoms forms 4 tt bands in Fig. 2b. 

Figure 2c schematically shows our experimental set- 
up. A back gate voltage Vg is applied across a thin layer 
of Si02 dielectric sandwiched by bilayer graphene and 
doped silicon to induce charge carriers in the sample. 
The gating efficiency is about 7.2 x 10^'^ cm^^/Volt [22]. 
For Stokes Raman scattering, ujstokes ~ i^l ~ ^s, where 
ujl and ws are frequencies of the incident and scattered 
light, respectively. The experiment is performed with the 
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FIG. 2; Bilayer graphene and the EFE-Raman setup. 

a, Bilayer graphene lattice with 4 carbon atoms At, Bt, Ai,, Bi, 
in its unit cell. Black arrows indicate carbon atom motion in 
the G band lattice vibration, b, tt bands of bilayer graphene 
near the K and K' corner points of the Brillouin zone. The 
low energy region is zoomed out to illustrate that the disper- 
sion is drastically different from single layer graphene [12]. c, 
Schematic drawing of the experimental setup. The two black 
lines in-between the gold contacts represent bilayer graphene. 
d. Optical image of the bilayer graphene device. The dashed 
blue circle is the position of the laser spot. 



sample mounted inside a variable temperature cryostat 
with optical access. Data shown in this work are taken 
at 12 K. 

The evolution of bilayer G band with the gate voltage 
Vg is displayed in Fig. 3. The EFE induced changes in 
the spectra are nearly symmetric about Vg = lOV (Fig. 
3a) . As in single layer graphene 0] , this symmetry, which 
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FIG. 3; Evolution of the G band of bilayer graphene with charge doping, a, G band at 5 different gate voltages. The 
arrows indicate the phonon peak positions. Spectra were taken at 12K in Helium gas environment, b, Energy and line-width 
of G band extracted from spectra in panel a. Two phonon anomalies are clearly resolved in the phonon energy evolution (blue 
squares) . 



determines the charge-neutral point of the sample, re- 
flects the underlying particle-hole symmetry in the band 
structure of bilayer graphene. Another important fea- 
ture in the data is that, the phonon bands have smaller 
line width at large charge doping with electrons or holes, 
indicating longer lifetime. This change of width is due 
to Landau damping of the G-phonon into electron-hole 
pairs when \Ef\ < tnjc/'i- at small charge doping. Sim- 
ilar effects were observed in EFE tuned Raman spectra 
in single layer graphene i^. 

While the observed G phonon line width evolution in 
bilayer is similar to that of a single layer, changes in the 
phonon energy are drastically different. In single layer 
graphene, G phonon frequency exhibits only one mini- 
mum as the EFE tuned charge density n passes through 
the charge-neutral Dirac point, and thephonon mono- 
tonically stiffens with increasing \n\ [1, 0, In con- 
trast, when charge carriers are added into the bilayer 
sample, fvjjQ first decreases for smaller |n| and then in- 
creases at larger doping (square symbols in Fig. 3b). 
As expected from the particle-hole symmetry in the sys- 
tem, two distinct minima are clearly resolved in hujcin). 
These minima indicate that G phonon stiffness changes 
from softening to hardening as carrier density increases. 

Non-monotonic changes of phonon energy as a function 
of carrier density is quite unusual as it has never been 
observed before. In chemically doped silicon 1^, 2^ and 
EFE doped monolayer graphene d, 0, [13], the phonon 
stiffness always changes monotonically with the carrier 
density. In this context, the double minima shown in 
Fig. 3b is quite anomalous. In the following, we ascribe 
that the observed minima are indeed manifestations of 



the general phonon anomaly described in equation (1) 
above. First, we note the similarities between Fig. Id 
and Fig. 3b, suggesting that the phonon energy min- 
ima correspond to Ep= ± hioojl in bilayer graphene. 
Furthermore, between the two minima, the phonon line 
width (red dots in Fig. 3b) is larger than outside of the 
minima, indicating that Landau damping of the phonon 
into resonant electron-hole pair transition (blue arrows in 
Fig. la,b) is allowed within this region Q. This is consis- 
tent with the fact that the anomaly positions correspond 
to the hole (left minimum) and the electron (right mini- 
mum) in the resonant electron-hole pair respectively. 

Since the phonon anomaly can only occur at a spe- 
cial Fermi energy, charge density difference between these 
two minimal points An^i « 6 x lO^^cm^^ is linked to 
the electronic band parameters that determine the low 
energy dispersion of the bilayer graphene. Within the 
tight-binding model [2^, we found from An^i the inter- 
layer Af, — Bi hopping energy 71 = 0.35 ± 0.06eV. This 
value is in reasonable agreement with 71 = 0.43± 0.03eV 
obtained from photoemission spectra in epitaxial bilayer 
graphene 24 1. 

We expect the phonon anomaly to show up in single 
layer graphene if the sample quality is improved. In pre- 
vious measurements [1, [ll| , single layer graphene samples 
contain large charge inhomogeneity, which yields electron 
and hole puddles of size = 3 — 10 x lO^^cm^^. This 
inhomogeneity corresponds to the mesoscopic Fermi en- 
ergy broadening bEp ~ 100 meV near the Dirac point in 
single layer graphene, a value large enough to wash out 
the anomalous phonon softening completely. In bilayer 
graphene, however, similar 8n results in much smaller 
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FIG. 4: Broadening of the phonon anomaly by charge density non-uniformity, a, Comparison of Ep ^ n relation in 
bilayer and monolayer graphene. The Fermi energy goes up much faster with charge density in the monolayer, b, Fits of the 
evolution of the G phonon energy with Ep. The best fit has charge density non-uniformity 5n — 1x10^'^ cm~'^. Blue squares 
are experimental data taken from Fig. 3b. c, Phonon anomaly A_Ba as a function of charge density non-uniformity in the 
bilayer (main panel) and the monolayer (inset). The anomaly is more robust in bilayer graphene. 



5Ef, because bilayer Ep changes much slower with n 
than the monolayer (n versus ^/n dependence in the low 
density regime), as shown in Fig. 4a. For this reason, the 
phonon anomaly is more robust and easier to observe in 
bilayer graphene. 

To be more quantitative, we evaluated numerically the 
effect of charge inhomogeneity on HuJciEF), assuming 
a Gaussian distribution of the charge density f{n) — 



"0 \ 

^ where no is the average charge density 
in the sample, 5n represents the size of the charge density 
non-uniformity. Figure 4b exemplifies several such evo- 
lutions for different 5n. As expected, the larger 5n is, the 
smaller is the phonon anomaly size AE'^ = Hujg{Ep—0) - 
huJG{Ep=hLUG/'^)- Figure 4c displays AEa as a function 
of 6n. Comparing with the experimental observation, 
we estimate that Sn ^ 1 x 10^^ cm~^ for 1 cm~^ AEa 
in a bilayer. This size of charge density non-uniformity 
agrees well with results obtained from other experimen- 
tal methods [25, 2^, 27 1. Note that similar size of charge 
density non-uniformity in single layer graphene is enough 
to smooth out the phonon anomaly completely (inset of 
Fig. 4c). 

In conclusion, we have observed the anomalous soft- 
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ening of the long-wavelength optical phonon in bilayer 
graphene. This striking effect is a manifestation of the 
predicted logarithmic divergence due to tunable resonant 
coupling of phonons with particle-hole pairs. The broad- 
ening of the phonon anomaly is attributed to large charge 
density non-uniformity in graphene, showing that EFE- 
Raman spectroscopy can access fundamental interaction 
effects near the charge neutral point of graphene layers 
even in the presence of an inhomogeneous charge density 
distribution. 



We thank I. Aleiner, D. Basko and A. Millis for helpful dis- 
cussions. We acknowledge financial support from NSF (CHE- 
0117752), the NYSTAR, and ONR (N000140610138). RK. 
acknowledges support from the FENA MARCO Center. A. 
R is supported by NSF (DMR-0352738) and DOE (DE-AI02- 
04ER46133). 



Note added. - After completion of this work, we noticed a 
theoretical paper [2^ on electron-phonon coupling in bilayer 
graphene, where it was confirmed that the long- wavelength 
optical-phonon energy again exhibits a logarithmic divergence 
as the charge density in the bilayer is continuously tuned. 
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